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ABSTRACT 

Context. One of the difficulties in extracting reliable information about the thermodynamical and magnetic properties 
of solar plasmas from spectropolarimetric observations is the presence of light dispersed inside the instruments, known 
as stray light. 

Aims. We aim to analyze quiet Sun observations after the spatial deconvolution of the data. We examine the validity 
of the deconvolution process with noisy data as we analyze the physical properties of quiet Sun magnetic elements. 
Methods. We used a regularization method that decouples the Stokes inversion from the deconvolution process, so that 
large maps can be quickly inverted without much additional computational burden. We applied the method on Hinode 
quiet Sun spectropolarimetric data. We examined the spatial and polarimetric properties of the deconvolved profiles, 
comparing them with the original data. After that, we inverted the Stokes profiles using the Stokes Inversion based 
on Response functions (SIR) code, which allow us to obtain the optical depth dependence of the atmospheric physical 
parameters. 

Results. The deconvolution process increases the contrast of continuum images and makes the magnetic structures 
sharper. The deconvolved Stokes / profiles reveal the presence of the Zeeman splitting while the Stokes V profiles 
significantly change their amplitude. The area and amplitude asymmetries of these profiles increase in absolute value 
after the deconvolution process. We inverted the original Stokes profiles from a magnetic element and found that the 
magnetic field intensity reproduces the overall behavior of theoretical magnetic flux tubes, that is, the magnetic field 
lines are vertical in the center of the structure and start to fan when we move far away from the center of the magnetic 
element. The magnetic field vector inferred from the deconvolved Stokes profiles also mimic a magnetic flux tube but 
in this case we found stronger field strengths and the gradients along the line-of-sight are larger for the magnetic field 
intensity and for its inclination. Moreover, the discontinuity between the magnetic and non magnetic environment in 
the flux tube gets sharper. 

Conclusions. The deconvolution process used in this paper reveals information that the smearing induced by the point 
spread function (PSF) of the telescope hides. Additionally, the deconvolution is done with a low computational load, 
making it appealing for its use on the analysis of large data sets. 

Key words, methods: data analysis, statistical — techniques: polarimetric, spectroscopic — Sun: magnetic fields, pho¬ 
tosphere 


1. Introduction 


Observations of the Sun from the Earth are always lim¬ 
ited by the presence of the atmosphere, which strongly dis¬ 
turbs the images. A solution to this problem is to place the 
telescopes in space satellites, which produce observations 
without any (or limited) atmospheric aberrations. Recent 
examples of the se atmosph eri c-free observations are the 


Hinode mission (Kpsjjg^gJjjl. 
tropolarimeter (SP, .Lites et a.l 


2007), especially the spec- 
120131) of the solar optical 
telescope (SOT.lTsuneta et~ 20081 ), and th e vector mag¬ 
netograph IMaX^Martfnez Pillet et al.ll201l|) on board the 
Sunrise balloon (Solanki et alJl2Q10l ). Although the images 
from space are not affected by atmospheric seeing, the op¬ 
tical properties of the instruments still limit the observa¬ 
tions. In the case of diffraction limited observations, the 
point spread function (PSF) establishes the maximum al¬ 
lowed spatial resolution, defined as the distance between 
two nearby structures that can be properly distinguished. 
In space observations, the central core of the PSF is typi¬ 


cally dominated by the Airy disk, which is a consequence 
of a physical limitation imposed by the diffraction. Even in 
a diffraction limited instrument, real PSFs have typically 
the shape of the Airy pattern, with very extended tails. 
These tails do not limit the spatial resolution but induce 
a dispersion of the light from different parts of the image, 
leading to what is commonly termed as stray light or dis¬ 
persed light. This effect produces that light observed in a 
spatial location at the focal plane is a combination of the 
light emitted in the object at relatively distant spatial lo¬ 
cations. Therefore, the contrast of the object (defined as 
the pixel-to-pixel variation of the illumination normalized 
to the average illumination) measured in the focal plane is 
typically smaller than the contrast in the original object. 

The presence of stray light is important both for imag¬ 
ing instruments and slit spectropolarimeters. A first suc¬ 
cessful attempt to correct for this effect i n imag ing instru¬ 
ments was carried out bv lMartinez Pilletl (|1992f) . where an 
analytical PSF with long tails was proposed and the im- 
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age was deconvolved from it following a least-squares ap¬ 
proach. Another method could be to consider the stray 
light contamination as the sum of two components: a spec¬ 
trally dispersed component and a parasitic component of 
the spectrally undisp ersed light caus ed by scattering inside 
the spectrograph dBeck et al.1 1201 111 . In addition, we can 
also find in the literature the mult i-object m ulti-frame blind 
deconvolution technique (jvan Noort et al.ll2005ll which has 
been used to correct for all the perturbing effects of the 
atmosphere and the instrument. The novelty of the latter 
method is that a very general functional form for the PSF 
is proposed and blindly estimated from the observations, 
together with the corrected images. 

The case of slit spectrographs is more complicated be¬ 
cause the images are not immediately available. Instead, 
an image is constructed by adding the different scanning 
steps of the slit at different times. Therefore, one has to 
make some assumptions about the stability of the object 
so that the reconstructed images can be used for decon¬ 
volution. Furthermore, given that the spectral resolution 
of slit spectropolarimeters is much larger than that of im¬ 
age instruments, the number of monochromatic images is 
much larger. This makes the application of any deconvo¬ 
lution scheme a much more computationally heavy task. 
For this reason, it has been customary to postpone the 
treatment of dispersed light to the inversion phase of the 
spectropolarimetric data. A physical model (atmospheric 
model+radiative transfer) is proposed to explain the ob¬ 
served Stokes profiles, and an ad-hoc contamination is lin¬ 
early added to account for the stray light. It is possible to 
find different ways of computing this contamination in the 
literature: from local approaches that compute the stray 
light_ in a box o f N x N pixels around the pixel of interest 
(lOrozco Suarez et al.112007 ,1 to global approaches that use 
an average Stokes / profile in the whole field-of-view. Global 
approa c hes are p referred t o loc_al ones for different reasons 
(lAsensio Ramos fc Manso Sainzlf201l i. essentially because 
the use of local approaches somehow make the inversion 
process uncontrollable. 

One way to proceed when carrying out an inversion of 
spectropolarimetric data is to simultaneously do the in¬ 
version and the image deconvolution. T he first e f fort in 
this direction has been carried out by Ivan Noortl ( 2012 1, 
in which a standar d inversion code for the Stokes param¬ 
eters (iFrutiger et al J l2000h is modified to simultaneously 
take the presence of the spatial coupling induced by the 
PSF into account. This represents the first realistic ap¬ 
proach to a full inversion of the Stokes profiles without an 
ad-h oc treatm e nt o f the stray light. The approach followed 
by I van Noorti ( 2012 1 is computationally complex. The rea¬ 
son is that the inversion of the Stokes profiles is carried 
out simultaneously with the spatial deconvolution, using a 
Levenberg-Marquardt algorithm, without a distinction be¬ 
tween the two processes. This algorithm needs to compute 
and invert a Hessian matrix th at is ver y l arge. Tty mini¬ 
mize the computational load. lRuiz Cobo fc Asensio Ramosl 
(l2013f) used a simplified approach in which the inversion is 
carried out in two steps: in the first step the spectropo- 
larimetric data is deconvolved from the known PSF using 
a regularization based on a Karhunen-Loeve transforma¬ 
tion or principal component analysis PCA (see, lLoevelll955h 
and then inverted using standard Stokes inversion codes. In 
this paper, we ex plain the techni que ( only briefly presented 
in lR.uiz Gobo fc Asensio Ramosl 12013 1 in detail, and use it 


to analyze spectropolarimetric observations of a quiet Sun 
magnetic element obtained with Hinode/SP. 

The magnetic field of the solar surface is structured in 
a wide range of scales. The largest are the sunspots and 
can reach sizes of Mm. If we move to smaller structures we 
find pores, plages, faculae, or small magnetic elements that 
could reach a size of 100 km. Early spectropolarimetric 
observations have revealed many of the fundamental prop¬ 
erties of these structures, finding, for example, that the 
magnetic field intensity of these magnetic el ements should 
be of the order of kG values (see, for instance. IStenflolll973i l. 
Combined with basic magnetohydrodynamical (MHD) 
theory, this deduct ion led to the development of the thin 
flux tube mo del (Steiner et al.l 119981: ISpruit fe Zweibell 
ll 9791: [Parker) 119761) as a fundamental element on the 
structure of the photosphere. As time passes, the de¬ 
signs of solar telescopes has improved, providing new 
observational data to analyze their apparent size, bright- 
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rouglily provide a complete picture of our current knowl¬ 
edge. All of these studies generally support the idea that 
virtually all of the small-scale structure in active and quiet 
network regions is composed of filamentary flux tubes of 
kG magnetic field strength. 

In the present work, we explain in detail the spatial 
decon volution techni que em ployed for the first time in 
[Ruiz Cobo fc Asensio Ramosl ( 201, 1). In addition, we also 
apply this method to quiet Sun Hinode/SP data for the 
first time, aiming to take advantage of the spatial decon¬ 
volution process to analyze the physical properties of quiet 
Sun magnetic elements. 


2. Observations and data analysis 

2.1. Observations 


The polarimetric da ta we used were a cquired with the spec- 
tropolarim eter (SP; Lites et aDl2013H on board the Hinode 
spacecraft dKosugi et al. 1120071! . We selected a data set with 
a field of view of 82" x 164" recorded at disk center on 
April 21 th , 2007 (see Figure [[]). The SP instrument mea¬ 
sures the Stokes vector of the Fe I line-pair at 630 nm 
with a spectral and spatial sampling of 2.15 pm pixel” 1 
and 0.16", respectively. The exposure time is 12.8 s per 
slit position, making it possible to achieve a noise level of 
7.0 x 10” 4 I c in Stokes V and 7.2 x 10” 4 I c in Stokes 
Q and U. Here I c refers to the mean continuum inten¬ 
sity in the granulation. This data set displays a signal-to- 


noise ratio ^12.8/4.8 times higher than that of the nor¬ 
mal Hinode/SP maps in which the exposure time is 4.8 
s. To calibrate the spectra, we averaged the intensity pro¬ 
file from the whole map and compared it with the F ourier 
transform spectromet er spectral atlas (|Kurucz et al.lll98~4l : 
[Brault fe Necke1ll987l) once it was convolved with the spec¬ 
tral PSF jyf Hinode. A_similar calibration was done by 
[Cabrera Solana et al.1 (l2007t ) (see, Eq. 1). We found a dif¬ 
ference between the intensity of both profiles which can 
be interpreted as parasitic light inside the instrument, de- 
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fined as “veil” in iRuiz Cobo fe Asensio Ramos (120131 . In 
the present observation, the estimated value for this veil, is 
<7=0.0357 referred to the continuum signal I c . The data is 
finally corrected as / /w (A) = (I or { A) - C)/(I c - C). We 
subtracted this value from the continuum intensity before 
normalization. Since we aim to analyze strong magnetic el¬ 
ements in the quiet Sun (see black and white structures in 
Figure [Q), we selected one isolated magnetic structure, far 
from the edges of the map (see the red square), with strong 
longitudinal field signals, to study its properties in detail. 

2.2. Deconvolution 

We face the problem of correcting two-dimensional spec¬ 
tropolarimetric data from the perturbation introduced by 
the PSF of the Hinode solar optical telescope. We obtained 
the two dimensional data by scanning a slit on the surface 
of the Sun and recording the information of the four Stokes 
profiles (/, Q, U, V) on each point along the slit for a set 
of discrete wavelength points around the 630 nm Fe I dou¬ 
blet. As a consequence, the data can be considered to be 
four three-dimensional cubes of images. We use the nota¬ 
tion 1(A), Q(A), U(A), and V(A) to refer to observed images 
at a certain wavelength A. In practice, given the scanning 
process, these are not strictly speaking images because each 
column of the image is taken at a different time. 

In general, in the standard image formation paradigm, 
the observed image I (for simplicity we focus on Stokes /, 
but the same expressions apply to any Stokes parameter 
given the linear character of the convolution operator) that 
one obtains in the detector after degradation by the atmo¬ 
sphere and the optical devices of the telescope at a given 
wavelength can be written as 

I=0*P + N, (1) 

where P is the PSF of the atmosphere+telescope in the im¬ 
age of interest, while O is the original unperturbed image 
that one would obtain with a perfect instrument without 
diffraction and without any atmospheric perturbation. The 
operator * is the standard convolution operator and the 
quantity N is the noise contribution in the image forma¬ 
tion produced at the camera. We assume that we are not 
in the low illumination regime and N follows a Gaussian 
distribution with zero mean and diagonal covariance ma¬ 
trix with equal variance crf^. The previous expression can 
be applied to individual monochromatic images, with po¬ 
tentially different PSFs for each wavelength. For simplic¬ 
ity, we make the assumption that the PSF is wavelength- 
independent, which turns out to be a very good approx¬ 
imation given the wavelength ranges that we are dealing 
with (less than 2.5 A in the Hinode/SP case) . Th e spe- 
cific PSF that we consider is described in Ivan N oortl (l2012il 
and o btained from the pupil specified bv ISuematsu et all 
(2008:1), which takes the entrance pupil of the telescope and 
the presence of a spider into account. Under the presence of 
uncertainties induced by the noise, any deconvolution must 
be treated under a statistical framework. Consequently, we 
only have access to the distribution of reconstructed im¬ 
ages. Using the Bayes theorem, the posterior distribution 
p(0|I, P) which describes the probability of the restored 
image given the observed image and information about the 
image-forming system is given by 


P(0|I, P) 


p(I|Q, P)p(O) 
p( I) 


(2) 


where p(I|0, P) is the likelihood or, in other words, the 
probability that an observed image I has been obtained 
given an original image O and the PSF. The quantity p( O), 
also named prior, encodes all the a-priori statistical in¬ 
formation we have for the original images (i.e., degree of 
smoothness, presence of large gradients, etc.). This prior 
contents the a-priori statistical information of the whole 
field of view of the image. This statistical information could 
be, for example, that the Stokes I is a function that is de¬ 
fined as strictly positive, or that the magnetic field has 
spatial correlation between pixels. Finally, p(I) is a nor¬ 
malization constant (termed evidence) that is independent 
of the unperturbed image. Under the assumption of un¬ 
correlated Gaussian noise in every pixel of the image, the 
likelihood can be written as 


N 


p )=n exp 


i 


(I k - (O * P) fe ) 2 


2 a N 


(3) 


where the product is done over all N pixels of the image, 
Ik represents the fc-th pixel of the observed image (in lex¬ 
icographic order), and (0*P)fc is the k-th pixel (in lexi¬ 
cographic order) of the original image convolved with the 
PSF. The previous formalism allows us to obtain the max¬ 
imum a-posteriori (MAP) image, i.e. the image that max¬ 
imizes the posterior distribution. In addition, we assume 
that the prior over restored images is flat (all images are 
equally probable), the MAP solution is equal to the max¬ 
imum likelihood solution. Assuming that the prior p( O) 
is flat is equivalent to not limit the deconvolution process 
with any a priori statistical restriction. This solution can be 
found by taking the derivative of the previous Gaussian like¬ 
lihood with respect to the original image O. The resulting 
equation can be solved iteratively with an algorithm known 
as the Gaussian version of the Richa rdson-Lucy (RL) algo¬ 
rithm dRichardsonlll972l : lLucvHl97lf 


O new = Oold + [I — O old * P] (g) P, (4) 

where the symbol (g> represents image correlation (which 
can be written in terms of the convolution operator). The 
previous iterative scheme does not guarantee positivity of 
the images even when the observed images are positive and 
thus somehow have to be forced for Stokes I. However, note 
that for the Stokes parameters, the pixel values can be posi¬ 
tive and negative. Additionally, since the RL deconvolution 
is a maximum-likelihood algorithm, it is sensitive to over 
reconstruction produced by the presence of noise. The most 
notable effect is the appearance of high frequency structures 
in the reconstructed image. To avoid this problem, it is cus¬ 
tomary to stop the iterative scheme before these artifacts 
appear. 


2.3. Regularization 

The most straightforward way to deconvolve two- 
dimensional spectropolarimetric data is to deconvolve the 
monochromatic images of the Stokes profil es i n a way s imi¬ 
lar to what is done with imaging data (e.g. Jvan Noort et al.l 
120051 1. This approach presents two drawbacks. First, the 
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Fig. 1. Continuum map, left, and Fe I 6302.5 A magnetogram, right. The presence of small and not very common magnetic 
patches indicates that this map corresponds to a very quiet Sun region. Red square marks the magnetic element studied 
in detail. 


number of deconvolutions one has to carry out is large. For 
instance, the spectropolarimetric data of Hinode SOT/SP 
contains 112 wavelength points. Second, many of these 
wavelengths contain practically no information in Stokes Q , 
U, and V. This is the case of the continuum wavelengths 
where, unless strong velocity fields are present, the polari- 
rnetric signal is expected to be zero. Therefore, one ends up 
in the difficult situation of having to deconvolve very noisy 
images. The nature of the RL algorithm then induces an 
exponential increase of the spatially high frequency noise, 
making the final images useless. 

In general, and as a consequence of the smoothing intro¬ 
duced by the PSF, some information is irremediably lost. 
This unavoidably transforms the deconvolution process into 
an ill-posed problem. In particular, a set of solutions with 
potentially diverging power in the high spatial frequencies 
are perfectly compatible with the observations. Standard 
spatial deconvolution techniques solve this dilemma with 
ad-hoc spatial filtering methods that avoid the divergence of 
high frequencies during the deconvolution process. Typical 
methods include setting a hard or soft threshold on the 
resulting modulation transfer function that avoids the ap¬ 
pearance of high frequencies in the resulting image. 

We pursue a regularized deconvolution. Contrary to the 
typical procedure in image deconvolution, the regulariza¬ 


tion that we propose acts on the spectral dimension and 
not on the spatial dimensions. We assume that the unper¬ 
turbed Stokes profiles at each pixel can be written as a lin¬ 
ear combination of the elements of a complete orthonormal 
basis formed by the eigenfunctions {(/>i(A)}. Consequently, 
any of the unperturbed Stokes profiles can be written as 

n x 

0(A) = 5>i&M, (5) 

i=1 

where N\ is the number of wavelength points along the 
spectral dimension. If only a few elements of the eigenbasis 
are enough to reproduce the unperturbed Stokes profiles, it 
is advisable to truncate the previous sum and only take the 
first N <C N\ eigenfunctions into account. Therefore, the 
unperturbed data is now described by a set of images u>i 
(that we term projected images), which are built by pro¬ 
jecting the Stokes profiles of each pixel on the orthonormal 
basis functions. 

Given that we have assumed that the monochromatic 
PSF is wavelength independent, the observed perturbed 
Stokes profiles are obtained after applying Eq. © 

Nx 

I(A)=^( Wi *P)<MA) + N, (6) 

i=1 
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Fig. 2. First eight eigenvectors obtained after the PC A decomposition of the observed Stokes parameters from Figure 
[TJ The corresponding eigenvectors for Stokes /, Q, U, and V are displayed from top to bottom and the order of the 
eigenvectors increases from left to right. 


where we have used the fact that the convolution operator 
only acts on the spatial dimensions. Because of the presence 
of noise, we can find the original unperturbed Stokes pro¬ 
files by computing the projection of the previous equation 
onto the orthonormal basis functions 

Nx 

(1(A), 0 fc (A)> = * P) (</>i(A), 0 fc (A)> + N, (7) 

i=i 

where {•, •} indicates the dot product of the two functions. 
The noise term still maintains the same statistical proper¬ 
ties because the basis is orthonormal, which allows us to 
simplify the previous expression leading to 

<I(A),0 fc (A))=w**P + N, (8) 

Consequently, the regularization process we used im¬ 
plies that we have to deconvolve the projected images (as¬ 
sociated to the basis functions <j>k{ A)) from the PSF and 
reconstruct the unperturbed image using Eq. (0. This de- 
convolution is done using the RL iteration of Eq. 0. 

The previous approach is valid for any set of orthonor¬ 
mal fu nctions tha t on e util izes to ex plain the Stokes profiles 
(e.g.. ldel Toro Iniesta fc Lopez Aristeil2003 fh However, the 
basis obtained after PCA is ideal in our case because the 
PC A decomposition transformation is defined so that the 


first principal component accounts for as much of the vari¬ 
ability in the data as possible, and each additional principal 
component in turn explains the largest variability in the 
data under the orthogonality constraint. Therefore, work¬ 
ing with PCA-projected images, we find that the real signal 
present in each pixel only appears associated with the first 
few elements of the basis set, while the remaining elements 
are used to explain the noise. Consequently, the influence 
of noise is largely minimized if we only focus on the maps of 
low-order coefficients. This is a huge advantage with respect 
to the wavelength-by-wavelength deconvolution. 

The procedure starts by building the M x N\ matrix 
of measurements, where the Stokes profiles (with the mean 
Stokes profile subtracted) are placed as the rows of a ma¬ 
trix for each one of the M observed pixels. This matrix (or 
equivalently its covariance matrix) is diagonalized using the 
singular value decomposition (e.g., iPress et all 1 19861) . The 
eigenvectors obtained after the diagonalization form a basis 
that is efficient in reproducing the observed Stokes profiles 
and only a few of them are needed. In principle, and accord¬ 
ing to Eq. 0, we should use the PCA eigenvectors obtained 
with the original Stokes profiles. Since we do not have ac¬ 
cess to those profiles, we have used the observed Stokes 
profiles to build this database. Unless the original profiles 
are radically different from those observed, we expect the 
eigenbasis to be efficient in reproducing the original Stokes 
profiles also. 
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Fig. 3. Comparison between the observed and deconvolved continuum maps, first column; and the observed and decon¬ 
volved Fe I 6302.5 A magnetograms, second column. The observed region corresponds to the red box in Figured] Finally, 
the four colored squares indicate the location of the Stokes profiles we examine in detail later. 


The first eight eigenvectors for the four Stokes profiles 
computed using all the pixels of Fig. |T] are shown in Fig. [2] 
For the quiet Sun, contrary to the case of an active region, 
the noise contribution appears in the first PCA eigenvectors 
(see the 7 th or 8 th eigenvector in the linear polarization pro¬ 
files). This property is a consequence of the predominance 
of low signal polarization profiles, that need to be carefully 
taken into account. For our analysis, we only selected the 
first eight families of eigenvectors for Stokes I and V, and 
the first four families of eigenvectors for Stokes Q and U. 


volution and the non linear inversion of the Stokes pro¬ 
files, the resulting code is computationally simpler. It has 
the adva ntage th at any of t he existin g inve rsion codes 
as SIR dRuiz Cobo fc del Toro Iniest.al li)92li. NICOLE 


_ l2014h Spinor rtFrutiger et all l2000h . 

or Helix-1- (lLagg et al.ll2004ll can be used directly. The only 
addition is a first step in which one has to carry out the 
spatial deconvolution with a code that we provide for free 
for the community in the web address indicated in the con¬ 
clusions. 


ISocas-Navarro et al 


2.4. Comparison with other approaches 

Iiithe case of the inversion strategy presented bv lvan NoortJ 
(|2012f) . the regularization is done through the selection of 
the number and position of nodes that describe the physical 
models. This has the advantage that the physical interpre¬ 
tation of the filtering process is easy: the method eliminates 
the high frequencies in the Stokes profiles that need per¬ 
turbations in the depth stratification of the physical mag¬ 
nitudes with more than 3 nodes in depth. Obviously, the 
number of nodes can be changed at will, but then the in¬ 
version of the whole map has to be repeated. 

Our filtering procedure consists, essentially, of a very 
similar indirect suppression of high frequencies in the im¬ 
age by a filtering of the spectral features. However, since 
we use an empirical complete basis set that, in princi¬ 
ple, reconstructs all the profiles in the field of view to the 
noise level, we do not eliminate any important spectral 
feature that is already present in the data. Additionally, 
given the effective separation between the spatial decon- 


3. Properties of the Stokes profiles 

We mainly focus in the analysis of the magnetic patch en¬ 
closure inside the red box of Fig. [T] This magnetic element 
displays the strongest polarization signals and it is one of 
the largest structures of the map. It shows an almost circu¬ 
lar shape of longitudinal fields that surrounds some gran¬ 
ules. 

3.1. Spatial properties 

The first property we study is the changes the deconvolu¬ 
tion process induced in the continuum and in the magne¬ 
togram images. We plotted, in the first row of Fig. [3] the 
original continuum map and the Fe I 6302.5 A magnetogram 
(calculated as the difference between the Stokes V profiles 
at ±86 mA from the rest center), while in the second row 
we show the same magnitudes after the deconvolution pro¬ 
cess. We see a clear enhancement in the contrast of the 
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Fig. 4. From top to bottom, each row corresponds to Stokes /, Q , U, and V, respectively. Inside each panel we plotted 
the original profile in black and the corresponding deconvolved profile in red. In addition, we marked each panel with a 
colored square to indicate its position on Figure [3] 


continuum features. The contrast, defined as the standard 
deviation of the brightness in the continuum normalized to 
its average value, increases from 7.6% in the original map 
to 11.9% in the deconvolved map. The other noticeable dif¬ 
ference is the increase in the sharpness of the structure, 
together with a slight decrease in its size. This effect is 
probably better observed in the magnetogram, where the 
structure becomes smaller in the deconvolved data, indicat¬ 
ing that the smearing produced by the stray light contam¬ 
ination affects the polarization Stokes profiles. In addition, 
the deconvolution process also reveals an enhancement in 
the Stokes V signals in the interiors of the magnetic el¬ 
ement. This enhancement indicates that the deconvolved 
Stokes V amplitude sometimes reaches a factor of 2 larger 
than its original amplitude. 

Finally, we observed that a reversal in the Stokes V 
polarity appears in the edges of some structures (e.g., see 
the small isolated black patch close to the yellow square at 
coordinates 3."5, 8."7 in Fig. [3]). Th ese ringing stru ctures 
may be similar to those reported by iBuehler et al .! (|2015tl 
where they found a magnetic field reversal related to mag¬ 
netic patches of opposite polarity and magnetic field inten¬ 
sities below 300 G. However, we believe that, in our case, 
the ringing effects might be generated by the deconvolution 
procedure itself. We are currently studying spatial regular¬ 


ization techniques (like total variation regularizations) to 
efficiently suppress these effects. 

3.2. Stokes profiles 

The Stokes profiles found inside the magnetic elements with 
strong longitudinal fields (see white and black patches in 
Fig. [3]) are almost identical in the different regions. The 
Stokes V profiles are nearly antisymmetric with high am¬ 
plitudes, between 6-15% (normalized to I c ), in the interior 
of each patch, and their amplitude decreases as we move 
toward the edge of the magnetic patch. Linear polarization 
signals, on the contrary, are negligible in these magnetic 
elements. Finally, Stokes I profiles display a slightly asym¬ 
metric red wing that indicates the presence of velocity gra¬ 
dients along the line of sight (LOS) in the center of the 
magnetic element. 

To show the effect of the deconvolution process in the 
Stokes profiles, we have selected four different pixels that 
display strong Stokes V signals. The position of these pix¬ 
els are indicated with colored squares in Fig. [3] The cor¬ 
responding original (black) and deconvolved (red) Stokes 
profiles are shown in Fig. [4] 

The Stokes / profile displays major changes, typically 
presenting a large increase (for granules) or decrease (for 
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Fig. 5. Comparison between the original and deconvolved area asymmetry, first column; and the original and deconvolved 
amplitude asymmetry, second column. The vertical solid line marks the position of the inverted pixels we analyze in 
following sections. 


intergranules) of the continuum signal with respect to the 
original signal. Most important, the deconvolved profiles 
show strong changes in the core of the line. In fact, we 
detect the effect of the magnetic field producing the split¬ 
ting of the er-components. We believe that this effect is real 
(and not an artifact of the deconvolution) and related to 
the presence of a magnetic field because in most cases the 
splitting is present in the Fe I 6302.5 A, which is the most 
sensitive line to the magnetic field, while it is barely visible 
in the Fe I 6301.5 A line, which is the least sensitive to the 
magnetic field. 

Concerning the linear polarization signals, they are al¬ 
ways below the noise level both in the original and in the de¬ 
convolved profiles. Circular polarization signals, on the con¬ 
trary, display strong Stokes V amplitudes, reaching more 
than 10% of I c in the original profiles. The deconvolution 
process affects the Stokes V profiles in, at least, two dif¬ 
ferent ways: it slightly changes the Stokes area asymme¬ 
tries and it abruptly changes their profile amplitude. As we 
mentioned before, we find cases where, depending on the 
surroundings of each pixel, the change of amplitude could 
reach up to twice of its original amplitude signal. 

3.3. Analysis of Stokes profiles asymmetries 

Area and sometimes amplitude asymmetries are related to 
the correlation between velocity and magneti c field gradi¬ 
ents along the line of sight (tilling et al. 1975 ). From the 
results of previous sections, it is clear that the Stokes pro¬ 
files change after the deconvolution process. For this rea¬ 
son, we intend to study how the Stokes V area and am¬ 


plitude asymmetries change. To analyze the Stokes V area 
and amplit ude asym met ries , we follow the definition used 
in lMartfnez Pillet et al.l (jl997l ). that is, the area asymmetry 
is obtained as 


SA = s 


E^(Ai) 

EiinAor 


(9) 


where the sum is extended along the wavelength axis and 
s is the sign of the Stokes V blue lobe (chosen as +1 if the 
blue lobe is positive and —1 if the blue lobe is negative). We 
choose the range of integration of the Stokes V signals from 
—0.43 A to 0.43 A around the Fe I 6302.5 A line center. 

Likewise, the amplitude asymmetry is defined as 


5a — — —(10) 
CLb T Ur 

where a?, and a r are the unsigned maximum value of the 
blue and red lobe of Stokes V. 

We calculate these quantities over the small fragment of 
the map contained the red square in Fig. [T| where we only 
considered pixels that show Stokes V amplitudes higher 
than 1% of I c . The results of this study are included in 
Fig. [5] The first row corresponds to the original data set 
while the second row shows the deconvolved data. The top 
left panel shows that the area asymmetry is positive for 
the whole magnetic structure, with maximum values close 
to 10%. This area asymmetry is higher at the edges of the 
structure and decreases to being roughly compatible with 
zero when we move to the center of the structure. This effect 
is also visible in the amplitude asymmetry of the original 
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data (top right panel), which presents values up to 20% 
at the edges. The sign of the amplitude asymmetry is also 
unaltered for the whole structure. 

The maximum value of the area and amplitude asym¬ 
metries do not appreciably change in the deconvolved 
data (bottom row). However, the magnetic element shows 
some pixels that have changed the sign of both asym¬ 
metries, mainly in the center of the structure and 
aro und some edges. Similar changes have been found in 
lAsensio Ramos et al.l (120121 1 where the authors concluded 
that the smearing produced by PSF induces changes in the 
Stokes V asymmetries which can be partially recovered us¬ 
ing reconstruction techniques based on the phase-diversity 
procedure. This change of sign, especially that of the area 
asymmetry, indicates that the deconvolution process allows 
us to detect a change in the gradient of velocity and mag¬ 
netic field along the line of sight. In the following section, we 
study how this information is interpreted by the inversion 
process to define some relevant properties of the magnetic 
concentration. 

4. Stokes profiles inversions 

To obtain physical information of the atmospheric parame¬ 
ters where the Fe I lines form, we carry out the inversion of 
the Stokes profiles using the SIR (St okes In versi on based on 
Response functions; iRuiz Cobo & del Toro Iniestal 1 1992 1 
code, which allows us to infer the optical-depth depen¬ 
dence of these atmospheric parameters at each pixel in¬ 
dependently. 

We analyze the magnetic element shown in Fig. [3] 
with two different approaches, which we describe in de¬ 
tail in the following. The main difference resides that, in 
one case, we set the microturbulence fixed and equal to 
zero while, in the other case, the microturbulence was set 
as a free parameter. The microturbulence correction has 
been applied in low spatial resolution observations to re- 
produc e the properties of the Stoke s I line core (for in¬ 
stance, l^festendor2|£laz^£t^I][2Q0ll) • Although we are us¬ 
ing Hinode/SP data, we were not sure if the spatial resolu¬ 
tion of these observations is high enough to avoid the use of 
the microturbulence as a free parameter in the inversions. 
Thus, we aimed to analyze the results of the inversions using 
these two configurations, that is, with and without, micro- 
turbulence. 

4.1. Solution 1. No microturbulence 

We employed a single magnetic component parameterized 
by seven nodes in temperature T(t 5 oo)E], five in the LOS 
component of the velocity vlos( t 50o)> five in the magnetic 
intensity B(tsoo), three for the inclination of the magnetic 
held with respect to the LOS 7 ( 7500 ), and one for the az¬ 
imuthal angle of the magnetic held in the plane perpendicu¬ 
lar to the LOS <^( 7500 )■ Variables such as micro- and macro¬ 
turbulence are hxed to zero and not inverted. At each iter¬ 
ation, the synthetic prohles are convol ved with the spectral 
transmission prohle of Hinode/SP dLites et al.ll2013 ). Since 
each node corresponds to a free parameter during the in¬ 
version, our model includes a total of 21 free parameters. 

1 The parameter 7500 refers to the optical depth evaluated at 

a wavelength where there are no spectral lines (continuum). In 
our case this wavelength is 500 nm. 


The number of nodes in B(t 50 o), 7(7500), and v L os(t50o) 
are necessary to reproduce the small area asymmetries 
ofjthe observed (see F ig. O circ ular p olarization prohles 
(lLandolh fc Landi DeglTnnocentilll996lf . but mainly to re¬ 
produce the complex shape displayed by the deconvolved 
Stokes / prohles. Given that the inferred physical param¬ 
eters could be reliant on the initial atmosphere, we min¬ 
imize this effect by inverting each individual pixel with 
100 different initial atmospheric models. These initial mod¬ 
els were constructed by randomly perturbing the temper¬ 
ature stratihcation of the Har vard-Smithsonian Ref erence 
Atmosphere (HSRA) model dGingerich et al.1 ll97lD . The 
rest of the physical parameters of the initial model (B(t 5 oo), 
7(7500) and v 7,05(7"500)) are extracted from uniform proba¬ 
bility distributions and considered to be independent of the 
optical depth. Out of the 100 solutions obtained for a given 
pixel, we keep the one that yields the best ht. 

4.2. Solution 2. Non zero microturbulence 

The second conhguration is similar to the first, but in this 
case we also set the microturbulence as a free parameter. 
We choose three nodes for the microturbulence, while the 
number of nodes for the rest of atmospheric parameters re¬ 
mains the same. Then, the total number of free parameters 
for this second configuration is 24. We also follow the same 
strategy of inverting each pixel with 100 random inversions 
and choose the one with the smallest x 1 2 value (see Eq. fill) 
from the different solutions. 

4.3. Examples 

Some examples of the results of the inversion of Stokes pro¬ 
files using these two configurations are plotted in Fig. [6] 
In this figure, we only focus on the deconvolved profiles, 
although we also inverted the original profiles to compare 
the results in the following sections. We can see that the 
resulting fits from the two different configurations (red and 
blue lines) are close to the deconvolved profiles (black lines). 
However, although the Stokes profiles are well fitted in gen¬ 
eral, the core of Stokes I is not well reproduced for the Fe 1 
6302.5 A line. We believe that a configuration with the pos¬ 
sibility to introduce in the atmospheric parameters abrupt 
changes in height is necessary to fit these profiles. None of 
the solutions can be discarded from the point of standard 
model comparison. However, the solutions using these two 
configurations have some physical contradictions that we 
will explain later. 

4.4. Vertical cut 

The results of the inversion of the pixels marked by the 
vertical solid line in Fig. [5] using the first configuration, 
with no microturbulence, are shown in Fig. [7] From left 
to right, we show the depth stratifications of the tempera¬ 
ture, LOS velocity, and magnetic field. The first row shows 
the atmospheric parameters retrieved form the inversion of 
the original profiles, while the bottom row shows the same 
results using the deconvolved Stokes profiles. Concerning 
the temperature, we find a strong enhancement between 
log 7500 = — 1 and log 7-500 = —2 inside the magnetized re¬ 
gion. This range is precisely the place of maximum sensitiv¬ 
ity of the pair of iron lines at 630 nm. Coincident with the 
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Fig. 6. Results from the inversion of the deconvolved Stokes profiles presented in Fig. [I] (see the red profiles). We plot 
the Stokes I profiles in the first row and the Stokes V profiles in the second row. We omit the linear polarization profiles 
because their signals are always below the noise level. We plot the deconvolved profiles in black, the results from the 
inversion using the first configuration in red, and the results from the second configuration in blue. We also show in each 
panel the difference between the deconvolved and the inverted profile using the same color code of the corresponding 
configuration, and a colored square indicating the position of the pixel on the map of Fig. [3] 


increase in temperature we also find strong magnetic fields. 
Additionally, the height of the atmosphere at which the 
magnetic field is strong increases when moving apart from 
the center of the magnetic structure. This can be taken as 
an indication that the magnetic field lines are starting to 
fan out. This opening of the magnetic lines is due to a de¬ 
crease of the gas pressure outside the magnetic element. 
Finally, this increase in the magnetic field and temperature 
is also consistent with a strong upward motion, while the 
plasma located outside the central part of the magnetic el¬ 
ement (at 0.7" and 1.6") shows downflow velocities. The 
values of the upflow velocities reach close to 4 km s -1 . 

In the case of the deconvolved data, second row of Fig. [7] 
the behavior of the physical parameters is slightly different. 
We find a temperature enhancement in the same region, 
although the spatial size of the structure is significantly re¬ 
duced. This is because the original structure is smeared by 
the spatial PSF and the deconvolution process has partially 
canceled this smearing. The magnetic field structure of the 
magnetic concentration is now different, showing a strong 
increase in the center of the structure. These higher inten¬ 
sity values are needed to reproduce the Zeeman splitting 
displayed by the deconvolved Stokes I profiles and occur at 
different heights at different pixels. In the central pixels, we 


find weak magnetic field values at middle heights (log T 500 
between —1 and — 2 ), and strong values at higher heights. 
Finally, the velocity of the deconvolved data shows a less 
smooth solution as compared with the original results. 

If we use the second configuration, that is, we set the 
microturbulence as a free parameter, we find the results 
shown in Fig. [ 8 ] The panel distribution is the same used in 
Fig.m but we added a column at the rightmost part of the 
figure that corresponds to the microturbulence results. 

Focusing on the original data (first row), the left most 
panel shows a minor enhancement of the temperature at 
the center of the magnetic element as compared with the 
first solution. The magnetic field shows essentially the same 
structure as Figure 0 concentrated field lines along the at¬ 
mosphere in the central part of the magnetic element, while 
these lines fan out when we move away from the center part 
of the structure. The third panel shows major differences 
between the results of the velocity along the LOS and those 
obtained with the first configuration. The LOS velocity in 
the central part of the magnetic element is close to zero 
and slightly moving downward, while for the rest of the 
region is slowly moving upward. The changes on the LOS 
velocity at the central part of the structure are responsible 
for the inversion of the area asymmetry sign found in the 
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Fig. 7. From left to right, temperature, magnetic field, and LOS velocity. The color code for the latter indicates in blue 
the upflowing material while the downflowing material is marked in red. The horizontal axis corresponds to the length 
of the vertical solid line in Figure [5] and the vertical axis corresponds to the optical depth. 


deconvolved data in some pixels, see Figure [5] Finally, the 
microturbulence is large in the central part of the magnetic 
element and decreases when we move away from the cen¬ 
ter of the magnetic structure. Its distribution resembles the 
magnetic field configuration. 

The second row of Fig. [ 8 ] shows the results of the inver¬ 
sion of the deconvolved profiles. As occurred with the first 
configuration, second row of Figure[Tl the spatial size of the 
structure is smaller in the deconvolved data case. Likewise, 
the magnetic field intensity is more intense in the inversion 
of the deconvolved data. In addition, we also find that some 
pixels display weak magnetic field values at middle heights 
(logr 5 oo = [— 1 , — 2 ]) and strong magnetic values at higher 
heights. Finally, the microturbulence also presents similar 
values between the original (first row) and the deconvolved 
data (second row), although with the same loss of spatial 
smoothness between the atmospheric stratification of adja¬ 
cent pixels as we found for the LOS velocity. 


4.5. The accuracy of the solutions 

It is clear from the previous discussion that two different so¬ 
lutions provide similar good fits. Given that both solutions 
give different configurations for some physical parameters, 
it is sensible to compare them using the \ 2 values to see if 
one of the two solutions is preferable. To do that, we use 
the following definition of the reduced % 2 : 


x = - 
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(ii) 


where, the superindex “obs” and “fit” designate the ob¬ 
served profile and the fitted profile. The quantity v is the 
degrees of freedom that corresponds to the difference be¬ 
tween the number of wavelength points and the number 
of free parameters. Finally, of is the noise variance of 
the original data and it has different values depending on 
each Stokes profile: 07 = 6.1 x 10 -3 , oq = 7.2 x 10 -4 , 
cry = 7.2 x 10 -4 , cry = 7.0 x 10 -4 . These values are normal¬ 
ized to the mean continuum signal obtained for the whole 
map. 

We show in Fig.[9]the ratio between the reduced y 2 val¬ 
ues of the solution with microturbulence over the y 2 values 
from the solution without microturbulence. We calculated 
this value for the inversions of the original and deconvolved 
data. We can see that they present similar values for the 
inversion of the original (black line) and the deconvolved 
(gray line) profiles as both lines are close to 1. In addition, 
the standard deviation of the ratio between the reduced 
X 2 values for all the inverted pixels is a = 0.4, and, con¬ 
sequently, both solutions are equivalent because the ratio 
values are mostly inside l±cr. We conclude then that the 
accuracy of both solutions is similar although the results 
for t he physical par a meters are, in many cases, different 
(see [Martinez Gonzalez et al.1 120061 . for a similar conclu¬ 
sion). However, we note that the first configuration is fa- 
vored w hen using the Bayesian Information Criterion (BIG; 
ISchwar zjl 1978T ). usually applied for approximate model com¬ 
parison. 
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Fig. 8. Same as Figured except we added the microturbulence results at the rightmost panel. 


Relciiitsjl y 2 



Fig. 9. Ratio between the x ' 2 values obtained from the 
second solution, with microturbulence, over the x ' 2 values 
of the first solution, without microturbulence. Black desig¬ 
nates this ratio obtained from the inversion of the original 
profiles, while gray depicts the ratio between both config¬ 
urations when we inverted the deconvolved profiles. These 
Stokes profiles correspond to the pixels marked with the 
vertical line on Figure [5j 

5. Discussion and conclusions 

The analysis strategy (decoupled deconvolution and in¬ 
version) we presented is very straightforward to use once 
the spatial PSF of the instrument is known. Since this is 
not the case for earth-based observations because of the 
presence of the fast-changing atmosphere, it is ideal for 
space-based observatories. We are convinced that the spa¬ 
tial deconvolution of 2D spectropolarimetric data prior to 


the i nversion o f the Stokes profiles (or during the inver¬ 
sion, Ivan Noortll2012D is the way to proceed in the near fu¬ 
ture. The stray-light contamination is correctly treated and 
o ne avoids t he da nger of incurring potential pitfalls (e.g., 
lAsensio Ramos fe Manso Saind 1201 111 . From the technical 
point of view, carrying out our analysis is straightforward 
with the standard tools freely available to any researcher. 
However, we make our IDL (and Python in the future) 
codcQ available online for everyone to use it, and also for 
the sake of reproducible research. 

The process of spatial deconvolution produces some 
changes in the Stokes profiles that we studied in detail. 
The first thing we found is an increase of the continuum 
contrast of quiet Sun regions from 7.6% for the original 
map to 11.9% in the deconvolved data, accompanied by a 
reduction of the size of the structures observed at contin¬ 
uum wavelengths and in the magnetogram. In addition, we 
detect a reversal in the Stokes V polarity at the edges of 
some of the analyzed magnetic structures. The validity of 
these polarity reversals has to be studied in future works 
because it can be an effect of over-reconstruction during 
the deconvolution process. 

The next noticeable difference that we found is the 
appearance of the Zeeman splitted er-components in the 
Stokes I profiles in regions of strong longitudinal field. The 
same pixels also display a sizable increase in the Stokes 
V amplitude. The analysis of the Stokes V area and am¬ 
plitude asymmetries revealed that the deconvolution pro¬ 
cess uncovers negative asymmetries in the central regions 
of magnetic structures. The absolute value of the maximum 
asymmetry does not appreciably change. 

We examined in detail a vertical cut crossing a mag¬ 
netic element. After the inversion of the Stokes profiles, we 

2 http://www.iac.es/proyectos/inversion/deconvolution 
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found good fits using two different inversion configurations. 
The first fit did not include the microturbulence as a free 
parameter. Using this configuration, we obtained enhanced 
temperatures on the core of the magnetized region, together 
with strong upflows. The magnetic field intensity displays, 
in the original and the deconvolved data, the schematic 
picture of a magnetic element whose magnetic lines are ex¬ 
panding with height. This s cena rio can be found in the liter¬ 
ature in theoretical works (iGrossm ann-Doerth et ah 2000(1 
as w ell as in observational results IlMartfnez Gonzalez et alJ 

l2012ll . 


The second configuration introduces the microturbu- 
lence as a free parameter and the LOS velocity distribu¬ 
tion completely changes; the plasma inside the magnetic 
element is slowly downflowing. We also found a tempera¬ 
ture enhancement at the locations of the magnetic element, 
albeit this enhancement is less important. The magnetic 
structure is similar to the one obtained with the first con¬ 
figuration. Finally, the microturbulence displays low values 
outside the magnetic structure although it could reaches up 
to 2 km s^ 1 inside the magnetic element. 

Concerning the differences between the inversions of the 
original and the deconvolved data, we found that the mag¬ 
netic structure becomes sharper and the smoothness be¬ 
tween consecutive pixels is less clear in the deconvolved 
data. The detected temperature enhancement increases in 
both cases while the magnetic field intensity strongly in¬ 
creases due to the necessity of reproducing the Zeeman 
splitting of the Stokes I profiles. 

From the two physical solutions obtained in the in¬ 
version of the Stokes profiles, we noticed that the ve¬ 
locity configuration of the second solution, see Fig. 0 is 
closer to the results of p revious studies rela t ed to the solar 
magn etic elements (see ISolanki et al.l 120061: Ide Wiin et alJ 
2009). However, the presence of enhanced microturbulence 
velocity inside the magnetic element is not expected in this 
part of the structure, which is roughly at rest. In addition, 
we expect the microturbulence velocity to be more impor¬ 
tant in outer parts of the magnetic element where we find 
the interface between the magnetic and non magnetic re¬ 
gions and where abrupt changes on the LOS velocity could 
be also found. 

It is also interesting to note the lack of coherence of the 
inferred LOS velocity in the results of the inversion of the 
deconvolved maps. This makes us think that the inference 
of enhanced velocities in the central parts of the magnetic 
structure is partially a consequence of the presence of the 
PSF. The inversion code fits the broadening of the Stokes 
I profile with this LOS velocity distribution. The deconvo¬ 
lution process removes this additional broadening and this 
compensation is not needed. 

In spite of using different inversion configurations, we 
could not perfectly fit the Stokes profiles. In particular, 
of interest is the Stokes I line core of the Fe 1 6302.5 A 
that displays the Zeeman splitting but with different in¬ 
tensities in each cr-component (see Fig. ©. If this is not 
an artifact of the deconvolution process, this indicates the 
existence of abrupt line of sight changes of the physical pa¬ 
rameters, probably with short height scales, at the line core 
formation region that could not affect to the Stokes V pro¬ 
files. These abrupt changes could produce a complex shape 
of the Stokes /, depending on the 7 r- and er-components 
shapes, while this effect would be not present in the Stokes 
V profiles. In fact, if the line is forming in this complex en¬ 


vironment, the presence of microturbulence on the second 
configuration of our inversion process is justified. 

Finally, we stress the fact that for the first time 
we have run the deconvol ution method presented in 
iRuiz Cobo fe Asensio Ramosl (|2013() on quiet Sun obser¬ 
vations. We have demonstrated that the strong magnetic 
elements display enough signal-to-noise ratio to reliably re¬ 
construct the information perturbed by the PSF and do 
not introduce many artifacts. We also point out that our 
approach allows us to pursue a trial-and-error study (un¬ 
avoidable when inverting weak signals) that is only possible 
because we have decoupled the spatial deconvolution and 
the inversion of the Stokes profiles. Although the true na¬ 
ture of quiet Sun magnetic elements is still far from being 
completely understood, we conclude that the spatial decon¬ 
volution of space solar observation will help to obtain more 
accurate results. 
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